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Abstract. Parton distributions in impact parameter space, which are obtained by Fourier transform- 
ing GPDs, are exhibit a significant deviation from axial symmetry when the target and/or quark is 
transversely polarized. From this deformation, we present an intuitive derivation of the Ji relation. 
In a scalar diquark model and in QED, we compare the Ji and Jaffe-Manohar decompositions of the 
nucleon spin. Using the MIT bag model, we estimate spectator effects through the presence of the 
gluon vector potential in the definitions of the quark orbital angular momentum. 

DISTRIBUTION OF QUARKS IN THE TRANSVERSE PLANE 

In the case of transversely polarized quarks and/or nucleons, parton distributions in 
impact parameter space [1] show a significant transverse deformation. In the case of 
unpolarized quarks in a nucleon polarized in the +x direction, this deformation is 
described by the _L gradient of the Fourier transform of the GPD E q [2] 

qq/p ,(xM)=J ^"* ±,A± ^° - A D-^/ ^-^^(^0 -Ai) 

(1) 

for quarks of flavor q. Since E q (x,0,t) also arises in the decomposition of the Pauli 
form factor F~ = /_j dxE q (x, 0, t) for quarks with flavor q (here it is always understood 
that charge factors have been taken out) w.r.t. x, this allows to relate the _L flavor dipole 
moment to the contribution from quarks with flavor q to the nucleon anomalous magnetic 
moment (here it is always understood that charge factors have been taken out) 

d* = J d 2 b ± q +x (x,b ± )b y = jjj^CO) = -^jXq/p- ( 2 ) 

Here e q K q / p is the contribution from flavor q to the anomalous magnetic moment of 
the proton. Neglecting the contribution from heavier quarks to the nucleon anomalous 
magnetic moment, one can use the proton and neutron anomalous magnetic moment 
to solve for the contributions from q = u,d, yielding K u / p ~ 1.67 and K d / p —2.03. 
The resulting deformation (\d q \ ~ 0.1 fm) of impact parameter dependent PDFs in the 
transverse direction (fig. 1) is rather significant and it is in opposite directions for u 
and d quarks. The sideways displacement of the center of momentum for each quark 
flavor from the origin provides a very intuitive derivation of the Ji-relation [3] for the 
contribution from quarks with flavor q to the nucleon angular momentum: Consider 
nucleon state that is an eigenstate under rotation about the i-axis (e.g. a wave packet 




FIGURE 1. Distribution of the j + density for u and d quarks in the _L plane (x — 0.3 is fixed) for a 
proton that is polarized in the x direction in the model from Ref. [2]. For other values of x the distortion 
looks similar. 



describing a nucleon polarized in the x direction with p = 0). For such a state, (T®°y) = 

= {T°y) and (T^z) = -(T*y), and therefore (T++y) = (T^z-T^y) = (J*). This 
result allows to relate the _L shift of the center of momentum for quark flavor q to the 
angular momentum P q carried by that quark flavor. The displacement of the _L center of 
momentum is a sum of two effects: 

• as discussed above, (T^ + y) for a quark relative to the center of momentum of a 
transversely polarized nucleon (Fig. 1) 

• however, already for a point-like transversely polarized spin 4 particle, the _L center 
of momentum (T ++ y) is shifted by i a Compton wavelength away from the origin 
(the center of the wave packet in the rest frame) 

In order to understand the 2 nd effect, i.e. the _L shift of the center of momentum for a 
_L polarized spin A particle, let us consider 'bag model' [5] type wave functions 1 for the 
wave packet of the target, but with a 'bag radius' R that will be sent to °° at the end 

Since XjT d z \y is even under y — > —y, iqy°d z q does not contribute to (T 0z y) = 
(iq (y°<9 z + f-d Q ) q). Using idoij/ = Ey, one thus finds 



■ J dWo z ay X f(r)(-i)dyf(r)y = ^- B J d\f{r). 



wy (4) 
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The precise shape of the wave packet does not matter as long as its size is sent to °o. 



In the limit when the bag radius goes to infinity E = M and / d 3 rf 2 (r) = 1 and therefore 
the 2 nd moment of _L flavor dipole moment for this wave packet reads 



T q ++ y) = (T°%) = I, (5) 



i.e. for a a spherically symmetric, delocalized (Dirac-)wave packet with J x = \ centered 
around the origin the _L center of momentum jfi{Tq + b y ) is not at origin, but at ^ ! 

This 'overall shift' of the nucleon _L center of momentum implies a contribution \(x q ) 
to (T ++ y) from quarks carrying momentum fraction \{x q ), which gives rise to the first 
term in the Ji relation. In addition, Eq. (1) implies a shift of the center of momentum 
for quark flavor q relative to that of the nucleon, by ^ / dxxE(x, 0, 0). Combining these 
two effects yields the Ji relation [3] 



J q = (T++b y ) = X -Jdxx [q(x) + E q (x,0,0)] 



(6) 



Since the famous EMC experiments revealed that only a small fraction of the nucleon 
spin is due to quark spins [4], there has been a great interest in 'solving the spin puzzle', 
i.e. in decomposing the nucleon spin into contributions from quark/gluon spin and orbital 
degrees of freedom. In this effort, the Ji decomposition [3] 

L L q q 

appears to be very useful, as not only the quark spin contributions Aq but also the 
quark total angular momenta J q = jAq + L z q (and by subtracting the spin piece also 
the the quark orbital angular momenta L z q ) entering this decomposition can be accessed 
experimentally through GPDs. The terms in (7) are defined as expectation values of the 
corresponding terms in the angular momentum tensor 



M ^' = ^^Z z q + Y f q 1 (rxiE>yq+ rx (e x fl) 



(8) 



in an appropriate nucleon wave packet. Here if) = id — gA is the gauge-covariant 
derivative. The main advantages of this decomposition are that each term can be ex- 
pressed as the expectation value of a manifestly gauge invariant local operator and that 
J z = jAq + L q can be related to GPDs (6) and is thus accessible in deeply virtual Comp- 
ton scattering and meson production and can also be calculated in lattice QCD. However, 
due to the presence of interactions through the vector potential in the gauge covariant 
derivative L z q does not have a parton interpretation. 

Recent lattice calculations of GPDs surprised in several ways [6]. First, the light 
quark orbital angular momentum (OAM) is consistent with L u — La, i.e. L u + Lj 0, 
which would imply that J g m j • 0.7 represents the largest piece in the nucleon spin 
decomposition. Secondly, L u —0.15 and L u m +0.15 in these calculations, i.e. the 
opposite signs from what one would expect from many quark models with relativistic 



effects, as we will also illustrate in the following section. While the inclusion of still- 
omitted disconnected diagrams may change the sum L u + L d , it does not affect the 
difference L u — L^. In Ref. [7], it was pointed out that evolution from a quark model 
scale of few hundred MeV to the lattice scale of few GeV might explain the difference. 

Jaffe and Manohar have proposed an alternative decomposition of the nucleon spin, 
which does have a partonic interpretation [8] 

\ = ^LA<7 + £^ + iAG + J^, (9) 

z z q q z 

and whose terms are defined as matrix elements of the corresponding terms in the +12 
component of the angular momentum tensor 

M+ 12 = ^ £^y 5 4+ +£?|(rx id) Z q + + e + ~ ij TrF +i A j + 2TrF +j (r x id) V. 

(10) 

The first and third term in (9),(10) are the 'intrinsic' contributions (no factor of rx) to 
the nucleon's angular momentum J z = +^ and have a physical interpretation as quark 
and gluon spin respectively, while the second and fourth term can be identified with 
the quark/gluon OAM. Here q + = jY Y + 1 i s the dynamical component of the quark 
field operators, and light-cone gauge A + = A + A z = is implied. The residual gauge 
invariance can be fixed by imposing anti -periodic boundary conditions A^(x^,oo~) = 
— A^(x^, — °o~) on the transverse components of the vector potential. 

Only the Aq are common to both decompositions. While for a nucleon at rest the 
difference in the Dirac structure between L z q and ££ z plays no role [9], the appearance of 
the gluon vector potential in the operator defining L z q implies that in general ££ z ^ L z . 
Other nucleon spin decompositions have been proposed in Refs. [10, 1 1]. 



QUARK OAM IN THE SCALAR DIQUARK MODEL 



In the scalar diquark model [12], the light-cone wave function for a 'nucleon' with spin 
t and the quark spin aligned and anti-aligned respectively reads [13] 



V/j, (x,k ± ) = (M + -W, 

+ 2 \ X / 



ki) 



\l/\(x,k ± ) 



k l + ik 2 



x 



0(x,ki) (11) 



with = ii \i +X 2 ■ H ere 8 is the Yukawa coupling and MlmlX are the masses 

of the 'nucleon'/'quark'/diquark respectively. Furthermore x is the momentum fraction 



carried by the quark and kj 



k±r represents the relative _L momentum. Using 



these light-cone wave functions it is straightforward to calculate a variety of observables 

t 

that appear in the context of nucleon spin physics. For example, since only i/r , carries 

~2 

(one positive) unit of OAM, the quark OAM according to JM is obtained as 



i dx i 



J 2 k_L 
16k 3 



(1-x) 



t 

1 

-7 



(12) 
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FIGURE 2. x distribution of the 'quark' OAM Jz^(x) (full) compared to L z q (x) from the unintegrated Ji 
relation (dotted) in the scalar diquark model for parameters A 2 = 10m 2 = 10X 2 . Both in units of ^§^2- 



where the factor (1 — x) is the fraction of OAM carried by the quark with momentum 
fraction x in a two-body system with one unit of OAM [9]. One can also use the above 
light-cone wave functions to calculate the GPDs that enter the Ji relation, to calculate 

1 /•! 

L q = - dx \xq(x) +xE(x,0,0)-Aq(x)} . (13) 

2 Jo 

Using a Lorentz invariant regulator, such as a Pauli-Villars regulator, one finds that the 
two definitions for quark OAM agree, i.e. L q = S£ q in the scalar diquark model, as was 
expected since L z q in the scalar diquark model does not contain a gauge field term. 

However, despite L q = J£ q , no such equality holds for the corresponding unintegrated 
quantities. If one defines L q (x) by Eq. (13) without the x integral [14], one does not 
obtain the OAM distribution that would be obtained from (12) without x integral (Fig. 2), 
which renders the interpretation of L q (x) as the quark OAM distribution questionable. 

ELECTRON OAM IN QED 

In QED, there are four polarization states in the ey Fock component. To lowest order, 
the respective Fock space amplitudes for a dressed electron with J z = +j read 

+ 2 + 1 X{i — X) +2~ X l—X 

W\ (x,k ± ) = (--mW,ki) Y t 1 ,(x,k 1 ) = (14) 

with 0(jc, k^_) = 4= kl — k2 2 . The label ±j represents the spin of the electron 

in the e 7 Fock component and ±1 that of the 7. 

Using these light-cone wave functions, it is again straightforward to calculate the 
OAM of the electron in the JM [8] decomposition. Including a Pauli-Villars subtraction 



with a heavy 'photon' with mass A one thus finds 



d\ ± 

167T3 



(1-x) 



+1-1 



vj/ 



x^o 



a 
An 



4, A 2 
-log — -- 

3 m z 9 



(15) 



Likewise, inserting the corresponding PDFs/GPDs from Ref. [13] into (13) yields 



T z 



1 1 

- dx [xq e (x) +xE e (x,0,0) - Aq e (x)] — 4 - — 

2 JO 47T 



4, A 2 7 
-log — + - 

3 m z 9 



(16) 



Both Jzf/ and L z e are negative, regardless of the value of A 2 (as long as A 2 > A 2 ). In 
the case of Sf z the physical reason is a preference for the emission of photons/gluons 
with the spin parallel (as compared to anti -parallel) to the original quark spin [15] — 
resulting more likely in a state with negative OAM. In fact, when x — > 0, this preference 
reflects the more general principle of helicity retention [16], which favors the lead parton 
(i.e. the parton carrying most of the momentum) to carry a spin as close as possible to 
that of the parent. It is also encoded in the evolution equations derived in Ref. [14]. The 
divergent parts of Jzf/ and L z e are the same so that their difference is UV finite 



^e Z ~L Z e 



A^oo 

An' 



Applying these results to a (massive) quark with J z = +j yields to 0{a s ) 



c/>z 



(17) 



(18) 



This result may have important phenomenological consequences. Recent lattice QCD 
calculations for GPDs yielded L z u < and L z d >0, which is opposite to the sign obtained 
in typical relativistic quark models — such as the MIT bag model. 

In QCD, the gluon spin is experimentally accessible, but the gluon OAM Jf z is 
not. On the other hand, the gluon (total) angular momentum J z appearing in the Ji 
decomposition is accessible, either indirectly (by subtraction, using quark GPDs), or 
directly, using gluon GPDs from lattice and/or deeply virtual J/xj/ production. Even 
though jAG and J z belong to two incommensurable decompositions of the nucleon 
spin, it is thus tempting to consider the difference between these two quantities, hoping 
to learn something about gluon OAM. Subtracting (9) from (7), it is straightforward to 
convince oneself that J z — jAG = ^g+Y^q {-^q ~ L z q ) , i.e. numerically J z — jAG differs 
from ££ z by the same amount that Y.q-^'q differs from Y* q L z q . In our QED example, with 
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d 2 k ± 
16k 3 
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(19) 



being the photon spin contribution, one thus finds (for A — > 0, A — > oo) 



5? z + — 



(20) 



As was the case in (17), ^ appears to be a small correction, but one needs to keep in 
mind that for an electron J z , Ay, and Jzfj are also only of order a. 

QUARK OAM IN THE BAG MODEL 

Since quark models have, to lowest order in a s , no vector potential, it makes perhaps 
more sense to identify the quark OAM from these models with ££ q rather than with the 
GPD-based L z . In Ref. [9] the difference A z e = Jz?/ — L z e was calculated to order a for 
a single electron in QED and the result then also applied to a single quark in QCD. 
However, in QCD quarks are never alone and the question arises regarding the effects 
from 'spectator currents' on the orbital angular momentum of each quark. 

In order to address this issue, we will in the following use the MIT bag model [5] to 
estimate @{a s ) corrections to the difference 

tf q = I} q -^ q = {j(rxgA) Z q). (21) 

The vector potential in (21) is calculated from the spectator currents, which are obtained 
by taking matrix elements in the corresponding ground state bag model wave functions. 
The vector potential resulting from these static currents is obtained by solving 

V 2 A a (r) = -?{r) = -^ g¥ }(f)a^ ¥ ,(r) (22) 

s' 

for each color component a and where the summation is over the spectators (here we 
pick the gauge V • A a = 0, but to 0{a s ) the result is actually gauge invariant — at least 
in the subclass of all gauges where all color components are treated (globally) SU(3)- 
symmetrically. I such gauges, matrix elements of operators of the type q^X a TqA a , where 
r is some Dirac matrix, and A a is calculated to &(oc s ), are proportional to the matrix 
elements of the corresponding abelian operators. It is thus sufficient to establish gauge 
invariance of q*r x Aq for abelian fields. The key observation is that the bag model wave 
functions contain no correlations between the positions of the quarks. Therefore, after 
eliminating the color in this calculation and introducing abelian currents, A5 factorizes 

into the density of the active quark \\f} (f)\\f s (f) times if x A) z = rA§. Writing the 
volume integral / J 3 r in cylindrical coordinates, one can isolate the only ^-dependent 
term r J 27r dtyA^ = § dr -Air) as a closed loop integral with fixed r and z. The closed 
loop integral is gauge invariant (its numerical value represents the color-magnetic flux 
through a circle with radius r) and so is the volume integral in which it enters. 
The contribution from a spectator with j z = s' to (q(r x A) z q) thus reads 

= ~\^t I rf3 ^ V ^^^^^(^)^^^(^)- (23) 

Note that A^ is independent of the angular momentum j z = s of the 'active quark', since 
y/j(r)i//i (?) = yr , (r)x\r_ \ if). However, it depends on the spin of the spectator since 

2 2 ~ 1 5 



the orientation of the vector potential entering (21) depends on the latter. For example, 
for s' = +j, one finds 

V/V)(?x af Wsl {r') = \^\ 2 Mkr f )Mkr f )2 X -^±^ (24) 
and hence independent of the bag radius [17] 

A^, =4*Mr/«ft7 M UUkr) +j \(kr)} 2 ^ ^ = T 0.78«, 

± 2 3 J r <R Jr'<R \r-r'\ f 

(25) 

When the active quark has s aligned with that of the proton, the two spectators must 
have opposite s' and their contribution to A cancels, i.e. A is nonzero only in those wave 
function components where the active quark has s = — \, in which case both spectators 
have s' = +j. As a result, A^ p is equal to twice A^_j times the probability to find that 

quark flavor with s = — \ (which is ^ for q = u and | for q = d, and hence 

A z , = -A z , = -0.052a s . A z = -A z , = -0.104a,. (26) 
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